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Abstract. The infrared properties of the high-temperature limit of Landau-gauge Yang-Mills theory are 
investigated. In a first step the high-temperature limit of the Dyson-Schwinger equations is taken. The 
resulting equations are identical to the Dyson-Schwinger equations of the dimensionally reduced theory, a 
three-dimensional Yang-Mills theory coupled to an effective adjoint Higgs field. These equations are solved 
analytically in the infrared and ultraviolet, and numerically for all Euclidean momenta. We find infrared 
enhancement for the Faddeev-Popov ghosts, infrared suppression for transverse gluons and a mass for the 
Higgs. These results imply long-range interactions and over- screening in the chromomagnetic sector of high 
temperature Yang-Mills theory while in the chromoelectric sector only screening is observed. 

PACS. 11-10. Kk Field theories in dimensions other than four - 11-10. Wx Finite-temperature field theory 
- 11.15.-q Gauge field theories - 12.38.-t Quantum chromodynamics - 12. 38. Aw General properties of 
QCD - 12.38.Lg Other nonperturbative calculations - 12.38.Mh Quark-gluon plasma - 14.70.Dj Gluons 
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1 Introduction 

The infrared structure of Quantum Chromo Dynamics 
(QCD) governs the dynamics of hadrons as well as the 
thermodymanic properties of hot and dense hadronic mat- 
ter and is thus of great current interest. Although our un- 
derstanding is far from being satisfactory much progress 
has been made recently using different genuinely non-per- 
turbative techniques. Such methods include lattice Monte- 
Carlo calculations, Dyson-Schwinger equations (DSEs), 
renormalization group methods, stochastic quantization, 
the use of topological arguments, and others. 

For QCD, being a gauge theory, it is expected that 
the description of strong interaction phenomena such as 
confinement contains gauge dependent aspects. Therefore 
investigations in a particular gauge will very likely not lead 
to a full understanding of non-perturbative dynamics but 
may nevertheless provide important information. It turns 
out that for the studies to be performed here the Landau 
gauge is advantageous [Tj due to its non-renormalization 
of the ghost-gluon vertex [2] . 

One of the most intriguing phenomena measured with 
extremely high precision is confinement, i.e. the absence 
of colored objects from the physical spectrum. This is a 
genuinely non-perturbative effect related to the infrared 
behavior of QCD. Monte-Carlo lattice calculations pro- 
vide clear evidence that the confining properties of QCD 
change above a critical temperature, see e.g. This has 
implications for relativistic heavy-ion collisions and the 
early stages of the universe. Since it is likely that confine- 



ment is generated in the Yang-Mills sector of QCD, the be- 
havior of a pure Yang-Mills theory may already reveal the 
qualitative mechanisms of the deconfinement transition. 
It is therefore especially interesting to study the proper- 
ties of Yang-Mills theory above the critical temperature 
as well as its temperature dependence. 

Here we will present the results of such a study. For 
technical reasons, this is currently done in the infinite- 
temperature limit, but investigations with finite-temper- 
ature corrections are ongoing. (The results of an earlier 
exploratory study are presented in T.) 

Complementary to this approach also a study at tem- 
peratures below the phase transition has been performed. 
The corresponding results will be presented elsewhere |S] . 
Both, these and the investigations presented here utilize 
DSEs inj, the equations of motion of the Yang-Mills the- 
ory. They therefore extend successful calculations in the 
vacuum of the Yang-Mills theory |3| and full QCD [H| to 
finite temperature. 

This paper is organized as follows: To make the pre- 
sentation self-contained we briefly review some aspects of 
confinement in Yang-Mills theories in section|21 In section 
Owe derive the DSEs to be solved and discuss the trunca- 
tions made. In section 0] we present analytic solutions for 
different momentum regimes while section exhibits the 
full numerical solutions as well as their comparison with 
lattice calculations. In sectional we discuss the thermody- 
namic potential and screening masses. Interpretations of 
the results with respect to different aspects are given in 
section 13 In section |H1 we conclude and also give an out- 
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look concerning ongoing activities. Technical details are 
deferred to five appendices. 



2 Confinement 

A major focus of the present work is the fate of confine- 
ment at temperatures far above the phase transition. It 
is therefore necessary to be able to extract information 
about confinement. To establish criteria for confinement 
we study the infrared behavior of the pertinent 2-point 
functions.^ Since calculations in Euclidian space-time are 
significantly simpler than in Minkowski space-time we use 
in the following three criteria which are applicable in the 
Euclidian casc.^ 

The first criterion is an empirical one. It allows to in- 
fer that a particle is confined irrespective of the dynam- 
ical origin of confinement. It is based on the fact that 
no Kallen-Lehmann representation of a particle exists if 
the corresponding propagator does not have a positive 
semi-definite spectral norm, i.e. for this "particle" the 
Osterwalder-Schrader axiom of reflection positivity 
is violated. It is then not part of the physical spectrum 
and thus confined, c.f. ref. QJI. To test this criterion we 
will investigate whether the Schwinger function A(t) is 
positive semi-definite, see section 16.21 below. In case the 
corresponding propagator D vanishes at zero momentum, 

lim D{q') = 0, (1) 

the respective Schwinger function cannot be positive, and 
thus positivity is violated. It is important to note that, 
since positivity violation is only a sufficient condition for 
confinement, particles may still be confined by other mea- 
ns. 

The two remaining criteria describe confinement mech- 
anisms. The Kugo-Ojima scenario puts forward the 
idea that all colored objects form BRST-quartets and thus 
do not belong to the physical state space. This scenario is 
based on a derivation requiring three preconditions which 
have (yet) to be proven. The first precondition is an unbro- 
ken BRST charge also for large and thus non-perturbative 
field configurations. The second is the failure of the clus- 
ter decomposition theorem. This requires that no mass 

^ In a SU{N) Yang-Mills theory 1-point functions are ex- 
pected to vanish even at non-vanishing temperatures due to 
the antisymmetry of the structure constants. No sign of a sym- 
metric color structure of vertex functions, although in principle 
possible for SU{N > 3) gauge theories, has yet been found 
In the following we therefore assume it not to be present. The 
dimensionally reduced theory, as the high-temperature limit of 
a four- dimensional Yang-Mills theory, then also only contains 
color antisymmetric vertex functions. 

Employing Euclidian space-time also facilitates compar- 
isons with lattice results. The transformation back to Min- 
kowski space-time is already non-trivial in the vacuum |1U| 
and will be addressed in section lS?^ For the case of equilibrium 
calculations this point is of minor concern since the statistical 
system is anyway described in Euclidian space. 



gap exists in the complete state space (although a mass 
gap has to exist in the physical subspace). The third re- 
quirement is an unbroken global color charge. Assuming 
the validity of these conjectures, the Kugo-Ojima confine- 
ment criterion is fulfilled in the Landau gauge if I14j 

lim q^Daiq^) ^ oo, (2) 
g^— >0 

where Dq is the propagator of the Faddeev-Popov ghosts. 
This scenario necessarily also implies the condition for 
the gluon propagator. 

In the Zwanziger-Gribov scenario entropy argu- 
ments are employed to show the dominance of field con- 
figurations close or on the Gribov horizon in field con- 
figuration space. It has two consequences that can be in- 
vestigated. This scenario predicts a behavior of the gluon 
propagator as in condition and of the ghost propagator 
as in condition In addition, it is argued that the gauge 
fixing term dominates the action. Therefore the ghost- 
loop-only truncation, introduced below, becomes exact in 
the infrared. Our results, as the ones obtained in the vac- 
uum 7,8 , support this picture in so far as and (0) are 
satisfied and that in the infrared limit the behavior of the 
gluon propagator is driven by the ghost loop. 

Also intuitively it is clear that a strongly divergent 
ghost propagator at zero momentum can mediate confine- 
ment. Such an infrared divergence relates to long-ranged 
spatial correlations. These are stronger than the ones in- 
duced by a Coulomb force since the divergence in momen- 
tum space is stronger than that of a massless particle. 



3 Dyson-Schwinger Equations 

3.1 Tensor structure of the gluon equation 

The DSEs pi6U17U18; form an infinite tower of coupled 
non-linear integral equations for the Green's function of 
a given theory. Thus it is necessary to truncate this sys- 
tem. The motivation for the specific truncation scheme 
used here will be given below. In the following we aim 
at a closed set of equations for the pertinent two-point 
functions. In Landau gauge and for the vacuum in 3-1-1 
dimensions these are the ghost propagator 

Dai,') - (3) 
and the gluon propagator 

D,,iq) ^ P,Aq)^ (4) 

where P^i.v{q) — ^iiv — <lyi,(lvlq^ is the transverse projector. 
Eqs. JSJ and Q define the dimensionless dressing func- 
tions G((7^) and Z{q'^), respectively. 

We will allow, but not require, ghost-loop dominance 
in the infrared. In the ultraviolet the one-loop terms are 
dominant as they are the leading terms in perturbation 
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Fig. 1. Graphical representation of the truncated system of 
vacuum DSEs. Wiggly lines represent gluons and dotted lines 
ghosts. Lines with a full dot are full propagators. Small dots 
denote bare vertices and open circles represent full vertices. 



theory. Thus we follow here ref. [7] and neglect one-parti- 
cle-irreducible two-loop diagrams. The derivation of these 
equations at zero temperature can be found e.g. in TTIITB]. 
The tadpole is not left as a free part of the equations but 
its non-perturbative behavior is fixed by the requirement 
that no divergences beyond those allowed by the Slavnov- 
Taylor identities should occur in the gluon equation. This 
reflects the behavior of the tadpole in perturbative calcu- 
lations. This truncation scheme is depicted in Fig. ^ 

Finite temperature is introduced in the DSEs using 
the Matsubara formalism ^|. This entails two indepen- 
dent tensor structures for the gluon, one of them being 
three-dimensional (3d) transverse and the other one 3d 
longitudinal 



(5) 



^ bo 



p- 




Fig. 2. Graphical representation of the truncated system of 
DSEs at non- vanishing temperature. Wiggly lines represent 3d 
transverse gluons, dashed lines 3d longitudinal gluons and dot- 
ted lines ghosts. Lines with a full dot are full propagators. Small 
dots denote bare vertices and open circles represent the full ver- 
tices which have to be constructed in the truncation discussed 
in the text. 

subsection 13.31 Pl and Pt are defined as [20] 



OQi, ^ 0^0001/(1 H — 5-j 



Pl^iv — Pfiv — PTfj.u- 



(7) 
(8) 



Graphically these equations are represented in Fig.|21 The 
solutions for small temperatures below the phase transi- 
tion are discussed in 



Therefore there are two independent dressing functions, 
Z and H.^ The variable q^ denotes the Matsubara fre- 
quency and q the spatial momentum. Note that Lorentz 
invariance, although not manifestly visible any more, is 
not lost in this formalism TT . The ghost propagator, as 
a scalar, does not acquire a second independent dressing 
function. It depends nevertheless also on go and q sepa- 
rately. Therefore three independent functions of two vari- 
ables each have to be determined. 

To obtain from the matrix equation for the gluon prop- 
agator two scalar equations for the dressing functions, it 
is contracted with two tensors ^ 

Pt,. = CPtm. + (1 - OiSf^u - (1 + ^)6M- (6) 

where ^ and C are two parameters. They are connected 
to the truncation scheme and will also be discussed in 



^ These are identically to Zt and Zl in respectively. The 
notation has been changed to better visualize their correspon- 
dence to the objects in the 3d theory. In |S] they are called Zm 
and Zo to emphasize their physical meaning. 



3.2 Dimensional Reduction 

In the next step we take the infinite temperature limit of 
the truncated DSEs. To this end we note that, neglecting 
any contributions from Matsubara frequencies different 
from zero, this yields an effective three-dimensional Yang- 
Mills theory with an additional adjoint Higgs. This Higgs 
field is the remnant of the Aq field in the four-dimensional 
theory. Thus the number of degrees of freedom is con- 
served. 

The structure of the projectors © simplifies in the 
case of vanishing Matsubara frequency, i.e. qo = 0. -Pf^^ 
becomes S^qSqi, independent of ^. It thus projects onto the 
time-time component of the propagator. Therefore, the 
3d longitudinal part relates to the Higgs and this sector 
corresponds to the chromoelectric sector of the original 
theory. The non-vanishing elements of P^^^ become 



p2 



(9) 



It projects onto the three-dimensional subspace. Therefore 
the gluons in the 3d reduced theory correspond to the 
chromomagnetic sector of the four-dimensional theory. 



4 



Axel Maas et al.: High- Temperature Limit of Landau-Gauge Yang-Mills Theory 



This dimensionally reduction amounts to integrating 
out the hard modes in tree-level approximation. Contri- 
butions from higher Matsubara frequencies can potentially 
generate additional terms in the Lagrangian of the di- 
mensionally reduced theory [22]. Lattice calculations and 
matching in the perturbative regime show that this is in- 
deed the case |22II2H| . Especially, a tree-level mass for the 
Higgs is generated. This also modifies the DSEs: the addi- 
tional terms occurring stem from the neglected Matsubara 
frequencies and once included they can be derived from 
first principles as follows. The Lagrangian which governs 
the 3d reduced theory is [TTII^ 

+ ^{DfcfJ'Dl^cj,^ + mlcjy^cj,^) + ^(/.'^ (10) 

with the field strength tensor F^^ and the covariant deriva- 
tive D"-}' defined as 

f;^ = d^^Al - d,Al - gsP'-AlAl (11) 
D';f ^d-'d^+g^r'-A';,. (12) 

A° is the gauge field, the ghost, c" the anti-ghost'', 
the Higgs field, g^, the dimensionful gauge coupling, 
mh is the Higgs mass, h is the Higgs self-coupling and 
j^abc ^jjg structure constants of the gauge group. A three- 
Higgs coupling is not present due to the antisymmetry 
of the structure constants and the remaining global color 
symmetry at tree-level. 

If appropriate, multiples of g^ are chosen as the fun- 
damental scale and all results in the employed trunca- 
tion scheme will be independent of Nc for gauge groups 
SU(iVc) in the following sense: In the chosen truncation 
scheme only the combination Ca(?| appears as reference 
to the gauge group, Ca — Nc being the second Casimir 
invariant of the adjoint representation of the group. Hence 
a change in Nc can always be absorbed by a change in g| 
as long as the ratio gl/ruh is kept fixed. 

The Higgs self-coupling h can be uniquely determined 
in the 3d reduced theory. Since the Higgs field is a com- 
ponent of the four-dimensional gluon field, no linear di- 
vergence of the Higgs self-energy can occur due to the 
Slavnov- Taylor identities of the four-dimensional theory. 
Introducing temperature, no novel divergencies arise liyv 
Implementation of this requirement fixes the Higgs self- 
coupling h in leading-order perturbation theory (see ap- 
pendix 0. 

All other constants occurring in the Lagrangian H10|l 
are effective constants which arise by integrating out the 
heavy modes. For any large, although finite, temperature 
T all remaining dimensionful quantities will scale as the 
temperature T. After rescaling, the precise dependence of 
33 on T is irrelevant for the solution of the DSEs and the 
issue will be postponed to subsection Ki.ll 

Note that the hermiticity assignment for the ghost is not 
the usual one but is still correct for Landau gauge |1I17| . 



The 3d reduced theory is not only superrenormalizable 
but finite, and all renormalization constants will be set to 
one in the following. 

The DSEs for the Yang-Mills sector are already known, 
see e.g. 17 . They are rederived in appendix IXI for com- 
pleteness. Since no tree level coupling between the Higgs 
and the ghost is present, the ghost equation will not be 
modified as compared to pure Yang-Mills theory. The gho- 
st and the gluon DSEs are given by 

xD^J{q)rf^''s^{^p,q,p-q) , (13) 

{Dh^Y\p) = S''\p^ + ml) + r^G;a6 ^ j.HH;ab 

xDf%q)r^bf{p + q,^p,^q) , (U) 
{D-Xl{p) = 5^\S,,p' - p,p,) + T^,^-^^' + T^J"'^" 

xD'a%P + Q)rf^''^'"'i-Q,P + Q,-p) 

+ ^/(0^"-""^^''^-^'-'^)^^'('^) 
xD^/,{p-q)r£-^^%^p,q,p~q) 

+ ^/(0^'''''^""^^'''-^'-'^)^"^'^) 
xD^fip-q)r^'f''''^f{-p,q,p^q) , 

(15) 

where are the tadpole contributions. The first label 
gives the equation where the tadpole contributes (G for 
gluon and H for Higgs) and the second the type of tadpole 
appearing. The index tl denotes the tree-level quantities 
which are explicitely given in eqs. H67|l - H72(l in appendix 
IXI This set of 3d truncated equations is diagrammatically 
displayed in Fig. |3| Note that all momenta in these equa- 
tions are three-momenta. 

It is important to note that the ghost-gluon vertex 
function r'^'^"*'°'"^(p, q, —p — q) becomes bare for vanishing 
ghost momentum p |21 , 

hm ^^'^^^(p, q, ^p~q)^ ig^P'^qu. (16) 

This identity especially ensures that we can choose a bare 
ghost-gluon vertex without altering qualitative results for 
the infrared behavior of the propagators. Indeed, at least 
in four dimensions, the qualitative nature of the infrared 
solution is independent of the detailed structure of the 
ghost-gluon vertex to a large extent [23] • Numerical con- 
sistency checks also find that the ghost-gluon-vertex does 
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Fig. 3. Graphical representation of the truncated system of 
DSEs in the infinite temperature limit. Wiggly lines represent 
gluons, dashed lines Higgs and dotted lines ghosts. Lines with 
a full dot are full propagators. Small dots denote bare vertices 
and open circles represent full vertices. 



not depart significantly from its tree-level form j25| . Thus 
we will use a bare ghost-gluon vertex in the following. To 
complete the description of our truncation we still have 
to choose the 3-gluon vertex and the Higgs-gluon vertex. 
As these choices are motivated by more technical consid- 
erations we will defer all details concerning these vertex 
functions to the appropriate sections describing the nu- 
merical solutions. 

To obtain equations for the dressing functions, the 
color indices are contracted assuming the tree-level color 
structure for the vertex functions. This assumption is sup- 
ported by lattice calculations ^ . The gluon DSE is trans- 
formed to a scalar equation by contracting it with the ten- 
sor ® . Dividing the DSEs by all trivial factors appearing 
on the l.h.s. these equations read 



47r2 



dqdeAT{p,q)G{q')Z{{p-q)^) 



(17) 



T 



HG ^ rpHH 



^ J dqde(N,{p,q)H{q^)Z{{p + qf) 

+N2{p,q)H{[p + qf)Z{q^)Y (18) 
GH I rpGG 



Z-\p^) = 1 + T^^ +T 



47r2 



dqde[R{p,q)G{q^)G{{p + qf 



+ML{p,q)H{q^)H{{jp + qf) 
+MT{p,q)Z{q^)Z{{p + qf) 



where 6 is the angle between p and q. For the truncation 
used in the present investigation, the integral kernels At, 
Ni, N2, R, Ml and Mt are given in appendix IbI^ 



3.3 Constraints on the solutions 

Due to the truncation of the DSEs, the Slavnov- Taylor 
identities are not fulfilled. This violation manifests itself 
in the appearance of spurious divergences in the gluon 
self-energy. These can be removed by adjusting either the 
tadpole terms or by projecting the gluon equation with 
the tensor ^ and choosing ^ = d = 3 in the equation 
for the 3d transverse gluon propagator [221 ■ The violation 
of gauge invariance does not manifest itself in the Higgs 
equation ffT^ of the 3d reduced theory. Thus this equation 
is independent of ^. Furthermore, the Higgs self-energy in 
3d is finite since the gluon propagator has only a logarith- 
mic divergence in four dimensions. Therefore the divergen- 
cies of the two tadpole terms in (|14|l have to cancel each 
other exactly. Nonetheless, as the Higgs has a finite mass, 
the tadpoles involving a Higgs propagator have a finite 
contribution which correspond to a finite mass renormal- 
ization. This has to be taken into account and will be 
discussed in detail in section^ 

Varying ( while selecting the tadpoles to cancel any 
spurious divergencies in the Yang-Mills and the Higgs sec- 
tor of the gluon equation (|15(l allows to explore different 
degrees of gauge violations. Especially the physical case 
of a transverse projected equation is of interest. If the re- 
sults do not depend qualitatively and only weakly quan- 
titatively on the value of it is justified to assume that 
the effect of the violation of gauge symmetry is small and 
under control, c.f. refs. |77) . 

There is another aspect of gauge invariance to be re- 
spected in non-perturbative calculations. As the Lorentz 
condition does not completely fix the gauge field, config- 
urations are over-counted. This is known as the Gribov 
problem |2H|- The region of gauge space which contains 
no copies and includes the origin is called the fundamental 
modular region, see e.g. There is, however, at present 
no local condition known which defines this region. It is 
contained in the first Gribov region including the origin 
and configurations for which the Faddeev-Popov deter- 
minant does not change sign. The restriction to the first 
Gribov region can be formulated at the level of dressing 
functions as boundary conditions to the solutions of the 
DSEs as 



G{q^) > 0, Z{q^) > 0, H{q^) > 0. 



(20) 



The condition on H does not follow from the Gribov con- 
dition in three dimensions but only due to the fact that 
it is a component of the gauge field in four dimensions. 
The condition 12U|1 does not completely solve the prob- 
lem since the first Gribov horizon encloses a larger part 
of gauge space than the fundamental modular region, and 



(19) 



^ For convenience, there have been slight redefinitions and 
changes of notation as compared to 
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therefore contains gauge copies [20] • Nevertheless, for the 
purpose of calculating propagators it is likely that this 
condition is sufficient for eliminating Gribov copies 15 . 
Note that, in lattice calculations, the Gribov ambiguity 
poses a more severe problem and makes it especially hard 
to extract the ghost propagator [ST], 



4 Infrared and ultraviolet behavior of the 
propagators 

4.1 Ultraviolet behavior 

The 3d reduced theory is also asymptotically free pS) . 
Thus the propagators and vertex functions reduce to their 
tree-level values for sufficiently large momenta. As the 3d 
reduced theory is finite, corrections are power-like. (In the 
case of the Higgs this is only true for k ^ rrih ) Since g| 
has the dimension of a mass and docs not enter in the 
loop integrals of eqs. ((T7|) . l(TH|) and ((T^ . already by di- 
mensional analysis, one infers that the one-loop contribu- 
tions are proportional to Thus all loop integrals are 
subleading in the ultraviolet compared to the tree-level 
contribution. An explicit verification of this behavior is 
given in appendix |0 Leaving aside the possibility of a 
finite renormalization, one has 

lim G{q^) = lim Z{q^) = lim H{q^) = 1. (21) 

— >oo — >oo — >oo 

For large enough the dressing functions are thus ap- 
proximately given by 



1 



(22) 



where the constants Ci turn out to be positive for all dress- 
ing functions. Hence all dressing functions approach the 
tree-level behavior from above. 



4.2 Infrared behavior 

In this subsection we will study the infrared behavior 
of the dressing functions employing the following general 
ansatze: 



G{q) = A,q-^\ 
Z{q) = A.q-^\ 



H{q) = A^q 



-21 



(23) 
(24) 
(25) 



For a pure Yang-Mills theory (without a Higgs) and a 
transversal projected gluon DSE this analysis has already 
been performed for three space-time dimensions in ref. 

m 

The assumption g = t = and Z = — 1, motivated by 
the naive idea of a free gas of quarks and gluons in the 
high-temperature limit, immediately leads to a contradic- 
tion due to divergencies in the ghost and the gluon DSEs. 



Thus no such 'Coulomb phase' can be realized when tak- 
ing into account condition Ijltill on the ghost-gluon vertex. 
Also g — and t = I = —1, corresponding to a purely per- 
turbative quark-gluon plasma with screening masses in the 
magnetic and the electric sector, immediately leads to a 
contradiction. This is also understood from the fact that, 
in order to obtain such a behavior, the tadpoles would 
have to dominate. However, for C = 3, the tadpole in the 
gluon equation does not contribute. Besides, it is hard to 
see how the vacuum corresponding to such propagators 
avoids the violation of Elitzur's theorem 34^ (for further 
details see ref. [SSI)- Hence the tadpole cannot be the in- 
frared leading term in the gluon equation and one must 
have i < — 1. Convergence of the integrals in the DSEs 
even requires t < —1, see appendix El This is in agree- 
ment with the reasoning [TFj for infrared dominance of the 
gauge fixing part of the Lagrangian (|10|l . The inequality 

t<-l (26) 

implies condition |^ and thus gluon confinement. On the 
one hand, this is no surprise as one expects confinement 
also in the three-dimensional Yang-Mills theory. On the 
other hand, it means that real-world chromomagnetic glu- 
ons are confined even at infinitely high temperatures! Note 
that the solution found in ref. also obeys the inequal- 
ity 

In appendix lElthe infrared limit of the DSEs is derived: 



A 



<'-^^/'-^-'lGT{9,t)AgA, 



m? + 5m? 

-7- = 1 + — ■ 

Ah y 



(27) 
(28) 
(29) 



where y — q^ and d denote the dimension. A subtraction 
in H27() has been performed and only the finite part is re- 
tained. The expressions for Iqt and Igg, originating from 
the ghost-self energy and the ghost-loop, can be found in 
appendix El In the Higgs equation, a finite renormaliza- 
tion of the mass has been allowed for. The mass renormal- 
ization is given in equation (|89|l and is discussed in |35j . 
Eq. H29(l possesses only one solution, namely / = — 1 and 



Ah = 



5m? 



(30) 



This indicates firstly a qualitative change in the high- 
temperature limit (as in four dimensions at T = one 
has t — I) and secondly the decoupling of the Higgs in 
the infrared. These observations agree with corresponding 
findings on the lattice |23) . 

The infrared ghost equation (|27() implies the relation 
(see also j^Sl) 



9 



-) 



The difference to the relation in four dimensions 



g =d^4 



-t/2 



(31) 



(32) 
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implies an additional power in momentum which exactly 
compensates the dimension of the effective coupling con- 
stant in three dimensions. A corresponding compensation 
is expected at high temperatures in four dimensions to 
cancel the straightforward temperature factor in front of 
the Matsubara sum. 

Due to the inequality (|26|l g is positive. Thus any solu- 
tion automatically satisfies Q and both, the Kugo-Ojima 
and the Zwanziger-Gribov conditions, are fulfilled. 

Using (jSU in Eq. (|2HIl one obtains 

- (2'^-'V^i^ + g){2 + die - 2) - 4g(C - 1) - 

,7T(d — 4(3) , , ^ 

X csc( ^ ) sm{ng)r{- + g) 



I Ghost exponent for different dimensions ~| 



'({d^lfgn 



l + d-2g 



1.(33) 



This equation has at least one solution for d > I and two 
solutions for d> 2, see Fig.0] Eq. I|33|) simplifies for d = 3 
to 

32g(l-5)(l-cot2(57r)) ^^^^ 



1 



(l + 2.g)(3-H2.g)(2-|-2.g(C-l)-C)' 

One of the two solutions in d = 3 is independent of the 
projection of the gluon DSE: 



(- --) 

4' 2' 



(35) 



while the other one depends on the parameter (. In the 
special case of a transverse projection, C = li one has 




Fig. 4. The upper panel shows the two solutions for the ghost 
exponent g for values of ( allowed by condition (1201 and inte- 
gral convergence, see Eq. (inSJ and ref. This excludes the 
second branch in d = 2. The dashed line is for d = 2, the solid 
line for d = 3 and the dashed-dotted line for d — 4. The lower 
panel displays the two solutions for the ghost exponent g as a 
function of d for C = 1. 



{g,t) = (0.3976,-1.2952) 



thus reproducing the results of ref. where only the 
case C = 1 has been treated. 

Condition H20|) . the inequality H26() as well as conver- 
gence of the integrals in the infrared restrict g in three 
dimensions to 

1<9<1. (37) 

Requiring furthermore a well-defined Fourier transform of 
the ghost propagator (at least in the sense of a distribu- 
tion) leads to the condition g < 1/2. This restricts the 
range of allowed C-values for the varying branch to 



(36) 

Determining these pre-factors turns out to be an essen- 
tial and unexpectedly complicated part of the numerical 
method. Thus Eq. (|39() is used to check whether a correct 
numerical solution is found (36| . 



1/4 < C < 3. 



At the lower boundary both solutions merge into one. 

Since in Eqs. (|27|l and (|28|l the prefactors Ag and 
only appear in the product A'^Az only this combination is 
determined by the infrared analysis: 



5 Numerical Results 

The numerical solution of the DSEs IHI), UHl) and (HHJ 
will be achieved in three steps: First, we study the ghost- 
loop only truncation where only diagrams with at least 
one ghost propagator are kept. Second, we include the 
gluon loop in the gluon equation. This truncation scheme 
(38) neglects all Higgs contributions and thus corresponds to a 
purely three-dimensional Yang-Mills theory. The last step 
is then to fully implement the complete system. A descrip- 
tion of the numerical method is given in refs. (23 



1 



AzA^ 



(47r)i 
24(s- 



1) 



(2 + 25(C-l)-C)/^(2-2.g)sin2(7rg) 



cos(27rg)(g-l).g2r(|-2g) 



5.1 Ghost-loop-only truncation 

Keeping only ghost loops and using a bare ghost-gluon 
vertex, the system of DSEs is completely specified. The 
results of the corresponding calculation with C = 3, for 
both the ghost and the gluon dressing functions, are shown 
(39) in Fig. In this case, only the 5 = 1/2 infrared solution 
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10"^ 



10 = 



k/g. 

Fig. 5. The solutions of the DSEs in the ghost-loop-only trun- 
cation scheme at = 3. The dashed line refers to the dressing 
functions of the gluon Z and the solid line to that of the ghost 
G. All quantities have been made dimensionless by dividing 
out appropriate powers of g\. 



exists. The gluon propagator exhibits a maximum located 
at kjgl « 0.25. 

The coupling constant is the only dimensionful quan- 
tity entering the DSEs. Thus all results can be represented 
through dimensionless variables when expressed in appro- 
priate powers of g^. Due to this fact, one can infer from 
Fig-Elthe results of the ghost-loop-only truncation for any 
positive value of the coupling constant. Even for a very 
small coupling constants strong non-perturbative effects 
for momenta smaller than g\ are clearly visible. Therefore, 
the infrared behavior of the pertinent 2-point functions is 
never perturbative. 

For ^ 3, a spurious divergence appears in the ghost 
loop of the gluon DSE. To cancel it, an appropriately ad- 
justed tadpole term is included. A subtlety required 
here comes from the fact that the subtraction due to the 
tadpole is only effective in the ultraviolet and not in the 
infrared. A detailed account of this subtraction is given in 
ref. |35j and the tadpole term T*^*^ is given in appendix 
IDl The solutions for the transverse projection {C, — 1) are 
displayed in Fig. El In this case, solutions for both sets of 
infrared exponents exist. As in the case C = 3, the solu- 
tions show no special features. 

To estimate the amount of gauge symmetry violation, 
the DSEs have been solved for different values of see 
Fig. 13 To this end we have varied C between and 4 for 
the set of half-integer infrared exponents. For the other 
branch, due to numerical uncertainties, C has only been 
varied from 0.25 to 2.55, the later corresponding Xo g — 
0.2708 and t = -1.042. Within the explored ranges the 
dependence of the solutions on ^, i.e. on the projection of 
the gluon DSE, is reasonably weak. 



5.2 Yang-Mills theory 

The gluon loop with a bare 3-gluon vertex in the gluon 
DSE contributes with opposite sign as the ghost-loop. This 
is correct in the ultraviolet and irrelevant in the infrared 
where it is subleading. At momenta of the order of g\ the 




Fig. 6. The ghost and gluon dressing function at = 1. The 
upper panel shows the ghost dressing function, the lower panel 
shows the gluon dressing function. The solid line denotes the 
solution for g = 1/2 and the dashed line gives the other solution 
branch with q ~ 0.4. 



gluon loop calculated with a bare 3-gluon vertex becomes 
dominant. This leads to zeros on the r.h.s. of the gluon 
DSE and thus to a violation of the Gribov condition l|20() . 
Therefore, within this truncation scheme, the tree-level 
3-gluon-vertex is not acceptable. Reasons might be that 
at momenta around (i) the bare 3-gluon-vertex over- 
estimates the true vertex, (ii) the bare ghost-gluon ver- 
tex underestimates the true one and/or (iii) the two-loop 
graphs are important. These possibilities will be explored 
in more detail in future work ,25, . However, as for the 
present study the specific reason for this shortcoming is 
of minor importance, we proceed by constructing an im- 
proved 3-gluon vertex. In this we are guided by earlier 
studies in four dimensions T'.'F' where minimal modifica- 
tions have been introduced to obtain the correct ultravio- 
let behavior. To respect the Bose symmetry of the 3-gluon 
vertex we model it by multiplying the bare 3-gluon vertex 
with an appropriate product of functions. This leads to 
the ansatz^ 



q,q + k, —k) 

-q,q + k,~k) {A{q')A{{q + kr)A{k')y , 
-- Z{q^)(G{q^)f^^l''3. (40) 



® Different ansatze, also ones violating Bose symmetry, have 
been employed, and it has been found that they do not lead to 
qualitative differences in the dressing function |85| . 
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Fig. 7. The ghost and gluon dressing function in ghost-loop- 
only approximation for different values of (. The upper panel 
shows the ghost dressing function, the lower panel the gluon 
dressing function. The solid lines represent the solutions with 
a ghost infrared exponent, g — 1/2, and the dashed lines the 
ones for the other set of infrared exponents. At the peak in 
the gluon, the middle lines correspond to solutions at — 1. 
The upper and lower line at mid-momenta give the solutions 
at = and 4 for the g = 1/2-branch and at = 0.25 and 
2.55 for the other branch. 



The function A is chosen to be a constant in the infrared. 
The additional parameter S does not change this behav- 
ior. On the other hand, it allows to smoothly interpolate 
between a large suppression and the tree-level value by 
tuning it from large positive values to zero. Stable solu- 
tions have been found for 6 > 0.114. In the following, we 
employ mostly 6 = 1/4 which yields a mild suppression 
of the 3-gluon vertex without effecting the stability of the 
numerical calculation. 

The resulting dressing functions, using C = 3, are sho- 
wn in comparison to the ghost-loop-only result in Fig. |S1 
For C 7^ 3 tadpole terms have again to be subtracted. 
Their construction is detailed in ref. j35| and the corre- 
sponding expressions are given in appendix ^ The re- 
sults for C = 1 ^iid (5=1/4 are shown in Fig. El The 
dependence on ( which provides a measure of gauge in- 
variance violation, can be inferred from Fig. 1101 and the 
one on S from Fig. 1111 The only qualitative difference in 
these functions, when compared to the ones of the ghost- 
loop-only truncation, is a maximum in the gluon dressing 
function which is now present for all employed values of 
The location and height of this maximum depends on the 
vertex construction. For most values of S this dependence 




Fig. 8. Solution of the Yang-Mills sector at = 3. The dashed 
curve gives the ghost-loop-only solution for comparison. The 
solid curve gives the full Yang- Mills solution at (5 = 1/4. Note 
the linear scale for Z in the lower panel. 



I Ghost dressing function 




Fig. 9. Solution of the Yang- Mills sector at = 1. The upper 
panel shows the ghost dressing function and the lower panel the 
gluon dressing function. The solid curve denotes the g = 1/2 
solution while the dashed curve displays the g ~ 0.4 solution. 
Both are at 5 = 1 /4. 
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Fig. 10. The ghost and gluon dressing of the pure Yang-Mills 
theory for different values of The upper panel shows the 
ghost dressing function, and the lower panel the gluon dressing 
function. The solid line gives the solution for g = 1/2 and the 
dashed line is for the other solution branch. At the peak in the 
gluon dressing function, the middle lines represent the solution 
a,t C, = 1. The upper and lower lines at mid- momenta give the 
solutions at = and C = 4 for the g = 1/2-branch and at 
C, = l/A and C, — 2.45 for the other solution branch. 



I Ghost dressing function" 
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Fig. 11. The ghost and gluon dressing functions of the pure 
Yang-Mills theory for different values of 5. The upper panel 
shows the ghost dressing function, and the lower panel the 
gluon dressing function. The solid line gives the solution for 
3 = 1/2 and the dashed line is the other solution branch. At the 
peak in the gluon dressing function, the middle lines represent 
the solution at 5 = 1/4. The upper and lower lines at mid- 
momenta give the solutions at (5 = 0.131 and 5 — 1 for the 
g = 1/2-branch and at 5 = 0.114 and 5 = 1 for the other 
solution branch. 



is weak, however. Nevertheless, we have to conclude that 
the behavior of the dressing functions is sensitive to tlie 
structure of the 3-gluon vertex and additional information 
about this vertex function is highly desirable and will be 
studied in the future. In the context of this paper we will 
rely on the comparison to lattice data (see subsection l5.4|l 
to justify our ansatz. 

Being confident that these technical issues do not in- 
validate our results we note that the infrared behavior 
of the gluon and ghost propagators in three-dimensional 
Yang-Mills theory is very similar to the ones within four- 
dimensional Yang-Mills theory. Differences are of minor, 
quantitative nature. This is not surprising since the three- 
dimensional Yang-Mills theory is expected to be also a 
strongly interacting, confining theory |32| . 



5.3 Including the Higgs 

Finally, we add the Higgs, i.e. we implement the full sys- 
tem of equations ((T7jl , (fTH|) and ifT^ . Although the Higgs 
loop in the gluon DSE provides a positive contribution it 
is not sufficient to allow for solutions with a bare 3-gluon 
vertex. Thus, the ansatz H4()|l will also be used for the full 



system. On the other hand, the tree-level gluon-Higgs ver- 
tex will be employed in this section. (In subsection l6 . 2l also 
a modified gluon-Higgs vertex will be studied.) 

As stated already, a tree level Higgs mass is induced 
when integrating out the higher Matsubara modes in the 
process of dimensional reduction '22^ . For the current anal- 
ysis, the origin and the exact value of the Higgs mass are 
not of direct importance. Hence this mass will be fixed 
to a value extracted from lattice calculations and we use 
fUh/gl = 0.8808 [231 in the following. 

As can be inferred from Fig.[5]three additional tadpole 
contributions arise when the Higgs is included, see also 
appendix^ Since the Higgs has a tree-level mass, two of 
the tadpoles can already have a non-vanishing finite part 
in leading-order perturbation theory, see appendixlO The 
self-energy of a Higgs field in a three-dimensional theory 
would be in general linearly divergent. This is not the case 
here. The Higgs, being a component of the gluon field in 
the four-dimensional theory, protects its self-energy from 
being divergent at the expense of fixing at leading-order 
perturbation theory the tree-level coupling h to be 

'. = -2.|j^. (41) 
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Fig. 12. The solution of the Yang-Mills system compared to 
the full system ai ( — 1. The solid line denotes the full so- 
lution at g = 1/2 while the dotted line is the corresponding 
Yang-Mills solution. The dashed line gives the full solution for 
the g ~ 0.4 solution while the dashed-dotted line is the corre- 
sponding Yang-Mills solution. 



see appendix [ni 

As can be seen from Fig. ^1 the influence of the Higgs 
on the Yang-Mills sector is small. This is in agreement 
with results from lattice calculations The resulting 
Higgs propagator, see Fig.^1 behaves similar to a massive 
tree-level propagator. 

Although the Higgs mass is fixed in the current set- 
ting to be rrih = 0.8808^3, the Higgs mass dependence of 
the gluon and ghost propagators, especially for decreasing 
values of the Higgs mass,^ would be of interest. However, 
already after a slight decrease in the Higgs mass, the so- 
lution ceases to exist. For masses below ruh ~ O.65I no 
solution is found (at least, for S < 1). 

The dependence of the Higgs on ( is very small, devia- 
tions between solutions for different values of C not exceed- 
ing more than a few percent. However, the Higgs dressing 
function, at momenta fc « 5I, is sensitive to the parame- 
ter 5, see Fig. Similar to what happens in the gluon 
dressing function, also the maximum of the Higgs dressing 
function increases with decreasing 5. This is understand- 
able from the fact that the Higgs self-energy depends on 
the strength of the gluon propagator. 



As expected, for large values of the tree-level Higgs mass, 
the solutions for the gluon and ghost propagators are indistin- 
guishable from the corresponding solutions of the pure Yang- 
Mills theory. 
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Fig. 13. The Higgs dressing function and propagator ai ( — 1. 
The solid line gives the solution for g — 1/2, the dashed line 
for g ~ 0.4, the dashed-dotted line denotes the leading-order 
perturbative result and the dotted line the tree-level behavior. 




Fig. 14. The Higgs propagator as a function of 5 at = 1. 
The solid line denotes the solution for g = 1/2 and the dashed 
line for g ~ 0.4. The higher the peak the smaller 5. The lowest 
peak corresponds to 5 = 1, the middle one to 5 = 1/4 and the 
largest one to S — 0.0862 and S = 0.061, respectively. 



Although the analysis reveals some dependence on the 
parameters C and 6 at intermediate momenta, it is im- 
portant to note that infrared properties are only weakly 
dependent on these superficial quantities. Only one of the 
two exponents depends mildly on ( as already discussed. 
In Fig.^lwe display the infrared coefflcient Ag as function 
of ^ or 5. (Note that Ah is fixed by the renormalized mass 
and Az depends uniquely on Ag.) Besides some depen- 
dence on S for small values of this parameter (where the 
solutions cease to exist, as well) the infrared coefficients 
are rather robust up to the edges of numerical stability. 
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Fig. 15. The dependence of the infrared coefficient Ag on ( in 
the upper panel and on S in the lower panel. Solid lines repre- 
sent the g = 1/2 solution and dashed lines the other solution 
branch. 

5.4 Comparison to lattice results 

Recently lattice results for the gluon and the Higgs prop- 
agator 122] have become available. As the gluon propaga- 
tor shows almost no sensitivity to the Higgs, not only in 
our but also in these lattice calculations, we additionally 
compare to the gluon propagator computed within three- 
dimensional lattice Yang-Mills theory |37| . 

As can be inferred from Fig. E| our results for the 
gluon propagator agree astonishingly well with the corre- 
sponding lattice results. As for the infrared behavior, the 
lattice results are in favor of the g « 0.4 solution. Note 
that, at momenta k ^ g^, our truncation scheme is not 
trustworthy and the better agreement of the 5 = 1/2 solu- 
tion with the lattice gluon dressing function is not conclu- 
sive. As displayed in Fig.Elthe DSE results for the Higgs 
propagator show clear deviations from the corresponding 
lattice results, the latter being closer to the leading-order 
perturbative results than to the DSE results. We will re- 
turn to a discussion of this point in subsection 16. 31 

6 Towards observables 

In a first step towards a calculation of observables from 
the propagators we determine an approximative form of 
the thermodynamic potential. Furthermore, we investigate 
screening masses via the calculation of the Schwinger func- 
tions from the propagators. In the following, the results 
forC = l,<5=l/4 and Ca ~ Nc — i will always be used. 
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Fig. 16. The gluon propagator and dressing function in pure 
Yang-Mills theory from the lattice and from subsection 15.2! 
The continuum-extrapolated values are from )23| . the others 
from |37l . The errors indicated are statistical. The solid line 
denotes the g = 1/2 solution and the dashed line the g ~ 0.4 
solution. Both are at = 1 and 5 — 1/4. 

6.1 Thermodynamic potential 

Knowledge of all Green's functions as functions of the tem- 
perature would allow to calculate the thermodynamic po- 
tential and therefore all thermodynamic quantities. Since 
not all of these are known, it is not possible to calculate an 
exact thermodynamic potential. Knowledge of the prop- 
agators is however sufficient to compute an approxima- 
tion to the thermodynamic potential, the Luttinger-Ward 
or Cornwall- Jackiw- To mboulis (LW/CJT) effective action 

This effective action can be extended to abelian and 
non-abelian gauge theories, see e.g. 39,40 . Its calculation 
is, however, prohibitively complicated when using con- 
structed instead of bare or exact vertices. Therefore we 
will only extract here some qualitative features using the 



Axel Maas et al.: High- Temperature Limit of Landau-Gauge Yang-Mills Theory 



13 



I Higgs dressing functional 




0.8 
0.6 
0.4 
0.2 





























p 






. . 1 , , . , 1 . . , . 1 . . . . 1 , . 


. , 1 . . . . 1 . , . . 




Fig. 17. The Higgs propagator and dressing function in the 
full theory from the lattice and from subsection l5.3l The lattice 
data are from |23| . The errors indicated are statistical. The 
solid lines give the g = 1/2 solution and the dashed lines the 
g « 0.4 solution, both at ^ = 1 and 5 — 1/4. The dotted lines 
represent the tree-level result while the dashed-dotted lines 
give the leading-order perturbative results. 



simplest approximation to the thermodynamic potential. 
It only depends on the propagators, and one has 

n ^ \d{G)T E / (Cf ( - HZ{q)) + {Z{q) 1) 
-fln(G(g))-(G(<z)-l) 

-^M|S) + Kf^-l))- (42) 

where d{G) = (5°° is the dimension of the gauge group. Hq 
is the tree-level Higgs dressing function 



(43) 



Working in the infinite-temperature limit we are eventu- 
ally interested in the energy density divided by T^. This 
motivates a rescaling of the integration momenta and one 
obtains 

^ = (44) 

where the dimensionless constant a depends only on the 
dimensionless ratio rrih/g'i ^^'^ thus becomes independent 
of temperature in the limit T — s- oo. In the simplest case, 
93 sKm)^ with the four dimensional coupling constant 
(74 depending on the renormalization scale /i. Herein enters 
the way in which the limit T — > oo is performed. When 
keeping 53 (/i) fixed by appropriately choosing /i(r), one 
has ~ J^QCD- Hence H44II scales as and does not 

contribute to the thermodynamic pressure significantly 
compared to the hard modes. The Stefan-Boltzman be- 
havior must then be obtained from the hard modes alone, 
requiring 



5|(m) 



(45) 



where b stems from the hard modes However, the soft 
modes may still contribute significantly to thermodynamic 
properties, especially the trace anomaly, near the phase 
transition '41'. 

The calculation of (|42|l turns out to be plagued by 
spurious divergences, since the DSEs with the modified 
3-gluon-vertex are no longer exact stationary solutions of 
()42|l . After carefully subtracting these a turns out 
to be within errors of 0(1), but significantly depending 
on the truncation and large numerical uncertainties. The 
best value for the full solution is 3.6 for the g ^ 1/2 solu- 
tion and 3.5 for the g « 0.4 solution, contributed mostly 
from the Higgs-sector. The difference of two phases can 
be extracted with better accuracy, as the leading spurious 
divergence cancels. The result is 



a 



g=l/2 



0.053 



(46) 



Thus the g ~ 0.4 solution is thermodynamically preferred, 
in agreement with the lattice results in subsection l5.4l 



6.2 Schwinger functions 

Screening masses are most directly extracted from the an- 
alytic structure of the propagators. To obtain access to 
these analytic properties we calculate the Schwinger func- 
tion related to the propagator D, defined as [10| 



Ait) 



1 



dpo cos{tpo)D{po), 



(47) 



i.e. the Fourier transform of the propagator with respect 
to (Euclidean) time. Note that this definition is indepen- 
dent of the dimensionality of the underlying theory. Neg- 
ative values for the Schwinger function can be traced to 
violations of positivity and therefore to absence of the par- 
ticle, represented by D, from the physical spectrum |10j. 
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Fig. 18. The Schwinger function of the gluon. The solid line 
gives the numerical result for g = 1/2 while the dashed line 
is for g « 0.4. The dotted and dash-dotted line denote their 
respective fits using the ansatz l|48|l . 



Solution 


A.ga 




/ 


5 = 1/2 


20.3 


64/27 


1.32511 




13.4 


64/27 


1.03148 



Table 1. The coefficients for the gluon fit lUSl . 



Using a sufficiently sophisticated FFT-algorithm to- 
gether with at least 512 or more frequencies^, it is pos- 
sible to calculate the Schwinger functions. The result for 
the gluon, presented in Fig. ^1 clearly exhibits positivity 
violations. This is in accordance with the Oehme-Zimmer- 
mann super-convergence relation and thus has been 
expected. Furthermore, the position of the zeros can be 
interpreted as the confinement scale. For the g — 1/2 so- 
lution, the zero occurs at g^t « 3.29 and at g^t w 3.98 for 
the other branch. This result is in agreement with recent 
lattice results |42| . 

To be able to perform an analytic continuation of the 
gluon propagator into the complex g^-plane wc first fit 
the Schwinger function. This is performed by using the 
ansatz^ 

zM- ,4S;,Ld i + TTkr,) (48) 

_ A.q^o+'{l+f+a^fq) 

for the dressing function. Az is the infrared coefficient de- 
termined previously, and a„ is the ultraviolet coefficient 
of leading-order resummed perturbation theory, as calcu- 
lated in appendix [CI The fit parameters for both solu- 
tions are given in tabled As demonstrated in Fig. llSlthe 
Schwinger function is fitted very well. The gluon propa- 
gator and dressing function are also reasonably well de- 
scribed by the fit, see Fig. 1191 

* Better are several ten thousands to a million. The results 
presented here have been obtained using roughly 5 ■ 10^ fre- 
quencies. 

^ In ref. |10| a parameterization with only a branch cut and 
no isolated pole provided a successful fit. Due to the different 
asymptotic behavior in four and three dimensions we use a 
different ansatz here. 



I Gluon dressing function ~j 




Fig. 19. Comparison of the fit-function (I48I I to the full so- 
lutions of the dressing function. The solid line indicates the 
g — 1/2 solution and dotted line its fit. The dashed line repre- 
sents the g ~ 0.4 solution and the dashed-dotted line its fit. 



Solution 


e 


fgt 





■m/g'i 


ff = l/2 


4.0199 


0.3545 


-0.67078 


1.4998 


subleading 


9.4493 


0.6736 


-0.18387 


2.561 


g^OA 


4.0697 


0.37793 


-0.63975 


1.5045 


subleading 


0.81188 






1.9223 



Table 2. The coefficients for the Higgs fit 1491 and the sub- 
leading one l|49|l and l|50|l . respectively. 



Similar to a meromorphic fit ansatz used in ref. 
we will describe the Higgs propagator and its Schwinger 
function, the latter being analytically calculated from the 
former, as follows 

Hfiq) = /^^•''r9l^ 4 , (49) 

A fit) = TrT-H^e-*"™"^'^) 
X (sin(0 + tm sin((/))) -|- sm{(f) — tm sin{(f>))) . 

As demonstrated in Figs. 1201 and 12 II these fits already de- 
scribe the Schwinger function quite well, but miss around 
20% of the propagator at zero momentum. This indicates 
that further massive modes are present. Indeed, for the 
5 = 1/2 solution, a further term of the form (|49|l has to 
be added. For the g « 0.4 solution, adding a term with 
one pole, 

HfiQ) = , (50) 



allows to improve the fit in this well. Both sublead- 

ing fits are not very accurate for large t, and the results 
have to be taken with care. They indicate, nevertheless, 
the existence of subleading contributions due to the pres- 
ence of further massive-particle-like contributions. The fit 
parameters of both solutions can be found in table [5] 

As can be seen from fig. El the result for the Higgs 
propagator deviates significantly from the lattice results, 
the self-energy being significantly overestimated. While 
the gluon Schwinger function is reasonably independent 
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I Higgs propagator^! 




I Higgs Schwinger function ^ 




I Higgs subleading Schwinger function | 




Fig. 20. The Higgs propagator in the top panel and its 
Schwinger function in the middle panel compared to their fits 
for g = 1/2. The solid lines represent the numerical solution, 
the dashed lines give the leading contribution and the dashed- 
dotted lines the first subleading contribution. The dotted lines 
underneath the solid lines give the sum of the leading and 
subleading contribution. The bottom panel shows the compar- 
ison of the numerical (solid) and the fitted (dashed) subleading 
Higgs contribution, see text. 
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Fig. 21. The Higgs propagator in the top panel and its 
Schwinger function in the middle panel compared to their fits 
for g ~ 0.4. The solid lines represent the numerical solution, 
the dashed lines give the leading contribution and the dashed- 
dotted lines the first subleading contribution. The dotted lines 
underneath the solid lines give the sum of the leading and 
subleading contribution. The bottom panel shows the compar- 
ison of the numerical (solid) and the fitted (dashed) subleading 
Higgs contribution, see text. 



of the truncation, this turns out not to be the case for the 
Higgs Sch"winger function. In order to sho"w this we study 
the sohition for a bare Higgs-gluon vertex suppressed via a 
scaling factor uj.^'^ First, we note that the gluon Sch"winger 
function, in contrast to the Higgs one, is not susceptible 
to such a change. 

Fig- mi displays corresponding results for different val- 
ues of Lo. The position of the first zero tends to increase for 
decreasing w. At a; w 1/4 any oscillation, at least "within 

In the Landau gauge, due to the transversality of the gluon, 
this is equivalent to modifying the tensor structure of the ver- 
tex. 



the available numerical precision, seems to be gone alto- 
gether. A fit to the Sch"winger function using reveals 
again additional structure at very small t "which cannot 
be captured by such a simple fit. As the fit also does miss 
some strength at zero momentum for the propagator this 
again indicates the presence of further massive contribu- 
tions in the propagator. 

6.3 Analytic properties 

Although the high-temperature limit of the four-dimen- 
sional Minko"wski theory is a genuinely Euclidian theory. 
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Higgs Schwinger function^ 




Higgs dressing function 
_ 1.4 




10"' 



Fig. 22. The top panel shows the Higgs propagator of the g ~ 
0.4 solution for various suppression factors cu of the Higgs-gluon 
vertex. The bottom panel shows the corresponding Schwinger 
functions. The lattice data are the same as in fig. 1171 The solid 
line is for uj = 1/4, the dashed line for u — 1/2, the dotted line 
is for w = 3/4 and the dashed-dotted line for a; = 1. 



it will be of interest for other applications to extract the 
analytic structure of the propagators investigated. 

The gluon propagator exhibits similar behavior for bo- 
th solutions but there are also some significant differences. 
The denominator of the ansatz (|48|) contains two factors, 
both of which could possibly give rise to a non-trivial an- 
alytic structure. The first part stems from the fit of the 
perturbative tail necessary to generate the maximum in 
the gluon dressing function. Since all fit parameters are 
positive, this factor does not give rise to a pole on the 
first Riemann sheet. However, it generates a pole on the 
second Riemann sheet, which will occur at l/a„/, that is, 
at Euclidian momenta. This pole does not have a physical 
interpretation, and may well be an artifact of the fit, since 
is tailored to generate the correct leading-order per- 
turbative behavior. Thus it generates most likely a struc- 
ture which has the Landau pole of perturbation theory 
on the second Riemann sheet. We expect that this pole 
vanishes when using a more sophisticated fit. 

The second factor generates a genuine isolated pole at 
(— (1 + f) / Az)~^^*' ■ This expression has only one value on 
the first Riemann sheet given by (—0.0933 -I- O.lGlSi)^! 
for the 5 = 1/2 and (-0.1462 + 0.1248i)g| for the other 
solution. In both cases a pole close to the origin is gen- 



Solution 


Order 


Pole/ffI 


ff = l/2 


Leading 


-0.5112 ± 2.191i 


ff = l/2 


Subleading 


-6.11972 ± 2.35777i 


g fa 0.4 


Leading 


-0.650078 ± 2.16822i 


5^0.4 


Subleading 


-3.69539 



" This is of the order of Aqcd when using the temperature- 
scaling of subsection l6.1l 



Table 3. Location of the poles of the Higgs propagator for 

UJ=l. 



erated with an imaginary part larger than the real part 
in one case. For the first solution, the pole is found for 
an angle significantly above 7r/4 while in the second case 
somewhat below. In addition, the first solution generates 
two more poles on two more Riemann sheets, while the 
second solution, with a (most likely) irrational exponent, 
generates an infinite number of further poles on an infinite 
number of Riemann sheets. In both cases, the residue is 
complex. In addition, there is a cut along the complete 
negative real axis starting at zero. In this way it is similar 
to the results in four dimensions |10) . 

From this analysis we infer that the gluon propagator is 
positivity violating, and thus satisfies the requirements of 
the Kugo-Ojima and Zwanziger-Gribov confinement sce- 
nario. 

On the other hand, the Higgs propagator very likely 
does not have a branch cut but a number of simple poles 
whose locations are given in table|21 The sensitivity to the 
Higgs-gluon vertex, however, necessitates further investi- 
gations before a firmer conclusion can be drawn. We only 
note here that for smaller values of the suppression factor 
w the Higgs propagator agrees better with lattice results, 
and its poles are then close to or on the real axis. 



7 Discussion 

Eqs. (fTTji . (fH^ and (fH^ describe the infinite-temperature 
limit of Landau-gauge Yang-Mills theory. The resulting 
gluon propagator exhibits the characteristic behavior of a 
confined particle. The Yang-Mills subsector is satisfying 
both, the Kugo-Ojima scenario ^ and the Zwanziger- 
Gribov scenario 101, thus describing a confined theory. 
These properties are very stable against the assumptions 
made: the result of confined gluons is nearly independent 
of the truncation and of the properties of the Higgs. Hence, 
in accordance with corresponding lattice results, the pres- 
ence of long-range chromoniagnetic forces is found. 

The chromoelectric sector is simpler in that it is close 
to a perturbative behavior. However, as the comparison 
to lattice calculations shows, at least next-to-leading-order 
perturbative effects or even genuinely non-perturbative ef- 
fects play a role. 

Combing these findings with the results from the vac- 
uum Uj and from the low-temperature calculations has 
impact on the understanding of the phase transition. The 
main difference between the low and the high-temperature 
phase is not primarily one between a strongly interacting 
and confining system and one with only quasi-free parti- 
cles. The chromoelectric gluons, whose infrared behavior 
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change from over-screening to screening, conies somewhat 
close to such a picture. The chromomagnetic gluons stay 
over-screened in the infrared and thus confined. The order 
parameter for the phase transition is then necessarily only 
a chromoelectric one. From the studies of Wilson loops 
it is known that only the temporal (electric) Wilson lines 
show a behavior typical of deconfinement while the spatial 
(magnetic) ones do not 03]. Note that the order param- 
eters used to study the deconfinement transition on the 
lattice are typically chromoelectric ones like the Polyakov 
lines |3] . Using the corresponding magnetic Polyakov lines 
we conjecture that almost no change will be found. This 
point of view is also supported by recent lattice calcula- 
tions 03] which observe even at 6Tc an over-screened mag- 
netic and a screened electric propagator, albeit on rather 
small lattices. 

The results presented here may have also consequences 
for the thermodynamic potential. In the analysis of the 
Higgs propagator indications for further particle- like poles 
have been found. These could contribute to the pressure. 
The way in which over-screened gluons contribute to the 
energy density but not to the pressure has been very re- 
cently discussed in ref. j41| . 

Another topic related to the investigations presented 
here is the recently established connection between four- 
dimensional Yang-Mills theory in Coulomb gauge and th- 
ree-dimensional Yang-Mills theory in Landau gauge 46 . 
The time-time component of the Coulomb-gauge gluon 
propagator is on the one hand directly linked to the static 
quark-quark potential and on the other hand to the con- 
nected parts of expectation values of the Faddeev-Popov 
operator. From this one concludes that the potential is ap- 
proximatively given by the ghost dressing function. The 
infrared behavior of the Coulomb-gauge ghost propagator 
in four dimensions and the Landau-gauge ghost propa- 
gator in three dimensions are determined from identical 
equations and one obtains identical infrared exponents. 
While the solution g ^ 0.4 generates a potential which 
behaves as ~ l/k^^^ and is thus a little less than linear 
with distance, the second branch generates a solution pro- 
portional to for small momenta, thus generating the 
behavior expected for a linear confining potential. 



8 Conclusions and Outlook 

We have analyzed the DSEs of three-dimensional Yang- 
Mills theory with and without an additional massive ad- 
joint Higgs field and solved them in a given truncation 
scheme. The investigated theory can be regarded as the 
high-temperature limit of a four-dimensional Yang-Mills 
theory. Besides the propagators we considered the cor- 
responding Schwinger functions. Although the resulting 
Higgs function behaves nearly perturbative it does not 
have a simple structure and experiences higher-order or 



even non-perturbative effects. The chromomagnetic glu- 
ons are over-screened. The Faddeev-Popov ghost propa- 
gator is infrared enhanced similar to the four-dimensional 
one. The corresponding long-range correlations thus im- 
ply confinement of chromomagnetic modes as they imply 
confinement of transverse gluons in four dimensions. 

We have given an approximate expression for the ther- 
modynamic potential. At the current stage it allows to 
discriminate which of the two sets of solutions is the pre- 
ferred one. Of course, it is our goal to extend our formalism 
such that physical observables such as the energy density, 
pressure and entropy can be calculated reliably. 

The next steps, however, also in connection with the 
results found in [S], will be to include the higher Matsub- 
ara frequencies to introduce the effects of a finite temper- 
ature into the system. This will hopefully allow us to de- 
termine the critical temperature and other aspects of the 
phase transition. Including quarks will allow to address 
the chiral phase transition. An extension to finite quark 
chemical potential is feasible 18 ,47], and thus investiga- 
tions of the QCD phase diagram based on a calculation 
of the infrared behavior of QCD Green's functions will 
become possible as well. 
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A Derivation of the Dyson-Schwinger 
equations 

The Dyson-Schwinger equations are the equations of motions 
for a quantum field theory. They can be derived jlRllTllj^ for 
the Euclidian version from 

((-^U^^+r(.))^[r(^)])^^^ = o (51) 

where (f) is the field variable, j the corresponding source term 
and a a generic multi-index. 5* is the Euclidian action and Z 
is the generating functional, 



-S[0°1+/ da:0''{a:)j»(a:) 



(52) 



that can be written as Z = exp(M^), W being a function of 
the external sources. The effective action, i.e. the Legendre 
transform of M^, is defined by 



Based on arguments |45| that fermion propagators will van- 
ish in the infinite temperature limit we do not expect that 
quarks can change anything in these considerations. 



w{f) = -nr) + J d''xf{x)r{x) (53) 

which implies 

- 5-7 . (54) 



18 



Axel Maas et al.: High- Temperature Limit of Landau-Gauge Yang-Mills Theory 



In the case of Grassmann fields, c, such as ghosts and fermions, 
two independent sources are necessary. This modifies the above 
to 

(^) ~ Srj'^(x) ^ (^) ~ ~ Sri"^x) 

= -r{u'',u'') + J d''x{fj''{x)c''{x) + c''{x)7j''{x)) (57) 

a / \ sr —a /\ sr 

where all derivatives with respect to Grassmann variables act 
in the direction of ordinary derivatives. 

The general procedure to obtain the corresponding Dyson- 
Schwinger equations is to calculate the expression 15111 for a 
given action and then perform once more a functional deriva- 
tive with respect to the field or with respect to the conjugate 
field in case of anti-commuting fields. The additional source 
term then yields the propagator while the right-hand-side of 
the equations are found by the derivative of the action. 

In general, propagators and their inverse are defined as 



S<P''{y)S4,'^(x) 
S4,>'(y)S4,'^(x) 



= D''\x-y)-^ 



(58) 
(59) 



while full n-point vertices are defined by a n-fold derivative of 
r . For the ghost-gluon vertex which is defined as 



&c'^{x)ic^{y)8Al,{z) 



rccA\ahc ( \ 



(60) 



the sequence of derivatives is relevant (in contrast to the gluon 
vertex functions). 

Using the identities 



5''W 



Sj^{x)Sri''(x) 

s^w s'-'r 



= - U'^zd-'w—r^ t . ^ (61) 

Sjl{z)5jf^(x) Sca{w)SAt(z) SriS(m)Sri'^(x) ^ ' 

^M4r,<:(a;) ~ I ^ (^'^t^ Sri" {x)Sc<' {z)^ Sri" (z) (^^2) 



S^W _ 

SAi(y)Sj^(x)Sf,(x) - 

-fdzdw 4^ ^ w (63) 

Sj^(x)Sjj,(z) SA'^(y)SAj,(z)5A-^{w) ^jt::{w)S]^{x)^ ^ ' 

it is straightforward (although tedious) to derive the DSEs in 
position space. Performing a Fourier transformation to mo- 
mentum space, where all momenta are defined as incoming 
and momentum conservation at the vertices is taken into ac- 
count, the DSEs in momentum space are derived. To make the 
appearance of the tree-level vertices explicit, it is necessary to 
add in general permuted versions of each diagram and using 
the (anti-)symmetry of the vertices in their respective fields. 
The ghost propagator equation reads 



^ r«;--^^''-(_g, p, q _ (p „ q)D'}? (g) 



the gluon one is given by 



(64) 



.rf-^:f'^\-q,p + q, -p) 



(27r)'' 



■r;^^t'%-P,q,P-q) 



A4>^;bef 



{~p,q,p-q) 



I 1 r d^q j^tl-.A'^ \abccl / „ „ „\ r^cd { \ 

Ti^^t^^^p \v-,-p,q,-q)D„p(si) 



H I T#Cfx"^'^(P. -p + q-k,k)- 
■Df^ip - q - k)Dll{q)DlUk)rC;t'^'{p -q~k, -p,q, k) 



,1 r d^_qd^ ptl\A'^ [acde j ^ „ „ ,^ ^ U\ Tidh I 
"•"2 J (27r)2<i n«70 



= (p,-g,g-fe-p,fc)D*(fc+p-g) 



Dli-\p) = 5-\5^.p^ -Pi^p.) 



■r%l,J-P,P + k, ~k)r^^^Jk + p-q,q,-k-p)- 
■D^(p + k)Dll(k)Df^{q) 

+h I ^ rf!^^""" ip, -p,i, ^q)D'H<i) 
--2 1 q-p-k,~q, k)Di\k) 

■DZip - q + fc)z?*(g)r/;^'^«'"^^(p - g + fc, ~p, q, -fc) 

+ I^n'^''''''-'-'''HP,q- P - k, -q, k)Dll{p + k-q)- 
■Df'(k)D^'' {k + p)I>* (g)r/'*''S^'=(p + k-q,q,-k-p)- 

■r,^^"--"''{-p,k + p,-k) , (65) 



and finally the one for the Higgs is 

D^\p)-^^S^\p^ + ml) 
+ I ■^r'J''^^'''^^^[^p - q,p, q)Dll{p + q) ■ 
■Df%q)r^'Hp + q, -p,-q) 
+ ^ I li^n'^^^'^'^'^'il, -q,P, -p)Dlt[q) 

-h I ^n'i''^'^''^'^%-P -1 + k, -k,p, q) ■ 
■D'^iF.{p + q-k)D%{k)D^\q)- 
■r^y"''"'\-k,p + q-k,-p,q) 
+ \ I k,p,-p + q- k)Dll ■ 

■D-^ip - q + k){q)D%{k)r^'^"'^^''{q,p- q + k,-p-k)- 
■D'"'{p + fc)r/'*''""^(-fc, ~p,p + k) 
+ 1 I ^rl'i'''^''-^'^'{-<l,1 + k,p, -p-k)Dll(q) ■ 
■D^\p + k)D'li{q + k)rf^-^''\k,q, ~q~k)- 
■D%{k)r^-^"'^^\-k,-p,p + k) 

H I (l^^*"^'^"''^'^, -p, 1, -q)D^^{q) 

-\! ^^r'^-'^'-'^^^'^ip, -p + q-k-k- -q)D'^^q) ■ 
■D^''{p -q + k)Df''(k)r'l''''-''^*'f[p -q + k,q, -p, -k) 
+ 1 I r"-'^' iP, -p-k + q,k, ~q)D^\p + k~q)- 

■D'^{q)Df''{k)rf-'''^'^^{q,p + k -q,-k~p)- 

■D^'^ip + k)r'*''''''^^{p + k, -p, -k) . (66) 
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Fig. 23. The DSEs of the propagators of the Yang-Mills theory 
with an adjoint Higgs given in I64L I65L and 16611 . The same 
conventions apply as in Fig.|3 



The tree-level vertices for the ghost-gluon-, 3-gluon-, 4-gluon-, 
2-gluon-Higgs-, 2-gluon-2-Higgs- and 4-Higgs- interactions ha- 
ve been used. They are given by 



II, 



(p,q,k) = igsf 



(67) 



mJp {P,<l,k) 



ritl-.A"^ '.abed / 7 j\ 

t^iwap (P,<l,k,l) 



+{p-q)pS^,.) 

j~\tl',A.(h'^ :abc / i \ ■ pabc/ j \ 

r,/ {p,q,k)=ig3f [q-k)^ 

2 c / p&ac rebd . read rebc\ 

ffaO^'A/ / +f f ) 



r 



tl^cf) \abcd 



{p, q, k, I) = 2h{SabScd + SacSbd + SadSbc), 



(68) 



(69) 
(70) 
(71) 
(72) 



where the momentum conserving (5-functions have boon sup- 
pressed. The complete graphical representation of 1164^ . H65|l . 
and 1661 1 is shown in Fig. 1231 



B Integral kernels 

The integral kernels in eqs. 1171 . 1181 and 11911 are obtained 
by using tree-level vertices and performing the corresponding 
contractions. The ghost kernel is given by 



At{k,q,e) = -- 



q^ sin^ie) 



(k^ + q'^ -2A:g cos 61)2 • 
The contributions in the Higgs equation are 

2<j2 sin^(6l) 



N^{k,q,9) = - 



N2{k,q,e) = - 



(fc2 +g2 ^2fcgCOSe»)2 

2sin3(6l) 



fe2 + g2 + 2kq cos e ' 



(73) 

(74) 
(75) 



The kernels in the gluon equations are finally 



R{k,q,e) = - 



((C-l)fcgcos(9)-g +Cg cos^(9)) sin 9 

ML{k,q,e) = 



((C-l)(fc +4fc(;cos9)-4g^+4g^Ccos^(fl)) sin 9 



4,kT(W+ip+2kq cos e) 



(76) 



(77) 



MT{k,q,e) 



4fc^(fc^+ij^+2fc<jcosfl) 

{k^ + 2q^)({(: - 9)fc^ - 4g2) + 8(C - 3)(fc2 + q^)kqcos 9 

+ (8(9* + (C + 7)k^ + 4(5C - l)fc'g') cos2(e) 
+4(4Cg2 + + 3)fc2) j,Qg3(g) _|_ 4^^2^2 cog4(g)\ ^^g^ 



The modified gluon vertex is introduced by multiplying Mt 

with iiinj. 



C Perturbative expressions 

Replacing all full quantities in the truncated DSEs l|17^ . H18|l 
and 11911 by their tree-level values, i.e. one for all dressing func- 
tions and the tree-level vertices H67I72|I instead of the full ver- 
tices, the standard perturbation theory to one-loop order is 
obtained. It is then possible to calculate the leading-order per- 
turbative dressing functions. For the ghost this becomes 



G{k)-^ = 1 



16fc 



(79) 



to leading order, independent on the presence of a Higgs. Note 
that a Landau pole at = 0{g^) is present. 

The calculation of the gluon self-energy is a little more 
tricky since, although finite in three dimensions due to a Slav- 
nov- Taylor identity, each single contribution is linearly diver- 
gent. Corresponding problems are most easily circumvented 
when contracting the gluon equation with the Brown-Penning- 
ton projector. Performing the calculation in pure Yang-Mills 
theory yields 

Z(fc)- = l-il|^. (80) 
Including the Higgs one obtains 



Z(fc)- = l-li||£A 



I 93CA I r, f 
^ 16-7rfc \ ^ k 



— p— ^csc 



1 + ' 



(81) 



Since non-perturbative effects already arise at order 1/k^ in 
the present truncation scheme, the only interesting part is for 
k ^ gi, ruh leading to 



Z(fc)-' = 1 



933 gA 

64fc ■ 



(82) 



Finally the Higgs self-energy is 

H{k)-' = 1 + ^ + ^- 



+ 



S3CA 1 



— 2k j arcsin( 



3) (83) 



where the third term is the tadpole contribution. The leading 
contribution is 

gICA 



ff(fc)-' = 1 



4fc 



(84) 
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The coupling constant h is determined by requiring that the 
sum of the tadpole kernels already generate a finite integral, 
i.e. independent of the regularization scheme. This is obtained 
by requiring I4H . 

Note that, in three dimensions, resummation has only ef- 
fects from order on since, by dimensional arguments alone, 
only tree-level expressions in the loops can contribute at or- 
der 33. Therefore H79^ . I|82^ and Ij84|l already constitute the 
resummed solution. This is confirmed by the numerical calcu- 
lations presented in section |3 

Note further that this also ensures that gauge symmetry 
is intact to one-loop order in the regime of applicability of 
leading-order resummed perturbation theory. The violation of 
gauge symmetry in this approach is, at least, not stronger than 
in ordinary leading-order perturbation theory. 



D Tadpoles 

The tadpoles used for the ghost-loop-only system read: 

yGG^ -f|^/rfade(i?c(fc,9)(G(g)G(fc + g) 

-Alq~^<^ik + q)~^^) + Ri[k,q){G{q)G{k + q) 

-Alq-^0{k + q)-''o (85) 

Here the kernel R has been split into its convergent, (^-in- 
dependent part Ro, its finite ^-dependent part iZs and its di- 
vergent part Rd as 

R^ Ro + Rs + Rd. (86) 
For the pure Yang-Mills theory this is altered to 

T^''^ -§^Jdqde(^RD{k,q)(G{q)G{k + q) 

~Alq~^^{k + q)-^^) 
+R3{k, q){G{q)G{k + q)- A^q'^^ {k + q)-^^ 

-1^ + MTD{k,q)Z{q)Z{k + q)j, (87) 

where Mtd contains the divergent part of Mr, including the 
alterations due to the dressed 3-gluon- vertex 1401 . For the full 
theory the additional tadpole in the gluon equation is given by 



T 



GH _ 9^Ca 



j dqde(^MLD{k,q)H{q)H{k + q)j . (88) 



(27r)2 

In the Higgs equation, the tadpoles are set to 

J.HG ^ J.HH ^ gfjCA m_ Sni^ 
fe2 471- ■ fc2 ' 



(89) 



A detailed account for the construction of the tadpoles is given 
in ref. |85| . 



the calculation of Igt and Iqg in H27|l and H28|l only. The lat- 
ter can be obtained straightforwardly when using the ansatz 
H25|l and the general formula 



J (27r)<i-^ ^ 



_ 1 r(-a-/3-j)r(f+a)r(j+/3) 4^^+^ , . 

^ '(I^ r(d-i-Q+^))r(-a)r(-/3) y v^^) 

valid for finite integrals, where x = q^ , z = {q~ kY and y = k^ . 
It yields 

Igg = ■ 

(47r)7 



2^g-^''(d-4g)(2-|-d(C-2)-4g(C-l)-C)r(d-2g)r(2g- j) 
(d_l)92p(l±i^)2p(j)2 

Note that the integral is convergent iff 
d~l d-2 

->g> — 



(91) 



(92) 



where the equality on the upper boundary requires the result 
to exist only in the sense of a distribution. 

Performing the same calculation for the ghost self-energy is 
more complicated. Since it is, in general, a divergent quantity 
Eq. H9()^ cannot be applied directly. It is necessary to regu- 
larize and then renormalize the expression. This can be done 
in a momentum subtraction scheme or via dimensional 
regularization. Here, the later has been performed by apply- 
ing the standard rules of dimensional regularization j48j . This 
then immediately gives a finite result. However, by doing so a 
divergent quantity has been removed which is formally elim- 
inated by setting —Z3 equal to this quantity. This procedure 
yields the same result as the momentum subtraction scheme. 
The range allowed for g in (I92I I renders Igt not only to have 
a divergence of logarithmic or linear order in even or uneven 
dimensions, but also quadratic or cubic divergences. Using a 
subtraction scheme, it would be necessary to include the next 
term in the Taylor expansion. On the other hand, dimensional 
renormalization directly yields 



Igt 

-29(4S(d_3)d+2l+29(l-|-9-d9))r(? 



93 Ca 
T 

(47r) 2 



-9)r(-9)r(29) 



(2-d+2g){d-|-29)r{f-29)r(9)r(f +g) 



(93) 



Note that this expression becomes negative already for values 
allowed by (I92L e.g. for g > 3/4 in three dimensions, thus 
reducing the allowed range and leading to the plots in Fig. |1] 



E Infrared expressions 



As already argued in section l4^ the only solution without fur- 
ther assumption is that of ghost dominance. This then requires 
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